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Abstract. We calculate the luminosity function of galaxies of the Early Data Release of the Sloan Digital Sky 
Survey (SDSS) and the Las Campanas Redshift Survey (LCRS). The luminosity function depends on redshift, 
density of the environment and is different for the Northern and Southern slice of SDSS. We use luminosity 
functions to derive the number and luminosity density fields of galaxies of the SDSS and LCRS surveys with a 
grid size of 1 ft -1 Mpc for flat cosmological models with Q, m — 0.3 and = 0.7. We investigate the properties 
of these density fields, their dependence on parameters of the luminosity function and selection effects. We find 
that the luminosity function depends on the distance and the density of the environment. The last dependence 
is strong: in high-density regions brightest galaxies are more luminous than in low-density regions by a factor up 
to 5 (1.7 magnitudes). 

Key words, cosmology: observations - cosmology: large-scale structure of the Universe 



1. Introduction 

The study of the distribution of matter on large scales 
is usually based on the distribution of individual galaxies 
or clusters of galaxies. An alternative is to use the den- 
sity field applying smoothing of galaxy or cluster distri- 
bution with a suitable kernel and smoothing length. This 
approach is customary in N-body simulations, where in 
each step a smoothed density field is evaluated. To the 
real cosmological data the smoothed density method has 
been applied in the study of th e topology of the galaxy dis- 
tribution by Gott et al. (|l986|). T he IRAS redshift survey 



was used by Saunders et al. ( 1991 ) to calculate the density 
fiel d up t o a distance 140 h^ 1 Mpc. Recently Basilakos et 
al. ( 2000 ) applied the same metho d usin g the PSCz-IRAS 
redshift survey by Saunders et al. ( [2000 ). In all three stud- 
ies a 3-dimensional spatial distribution was found. Due to 
the small volume density of galaxies with known redshifts 
a rather large smoothing length was used. This was suf- 
ficient to investigate topological properties of the galaxy 
distribution in the first case, and to detect superclusters 
of galaxies and voids in other cases. 

In this paper, we shall calculate the number and 
luminosity density fields based upon the Early Data 
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Release (EDR) of the S loan Digital Sky Survey (SDSS) by 
Stoughton et al. (2002) and the Las Campanas Redshift 
Survey (LCRS) by Shectman et al. ( |l996| ) . The SDSS 
Early Data Release consists of two slices of about 2.5 
degrees thickness and 65 — 90 degrees width, centred on 
celestial equator, LCRS consists of six slices of 1.5 de- 
grees thickness and about 80 degrees width. The number 
of galaxies observed per slice (over 10.000 in the SDSS 
slices and about 4,000 in LCRS slices) and their depth 
(almost 600 ft, -1 Mpc) are sufficient to calculate the 2- 
dimensional density fields with a high resolution (here h 
is the Hubble constant in units of 100 km s _1 Mpc -1 ). 
Using high-resolution number or luminosity density maps 
with smoothing scale of the order of 1 hr 1 Mpc it is pos- 
sible to find density enhancements in the field, which cor- 
respond to groups and clusters of galaxies. Using a larger 
smoothing length we can extract superclusters of galax- 
ies as done by Basilakos et al. ( 2000| ). In calculating the 
density fields we can take into account most of the known 
selection effects, thus we hope that the density field ap- 
proach gives additional information on the structure of 
the Universe on large scales. Clusters of galaxies from the 
SDSS were extracted previously by Goto et al. (2002)] us- 
ing the cut and enhance method. Loose groups of galaxies 
from LCRS were found by Tucker et al. ( J200C| ) . In accom- 
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panying papers by Einasto et al. (2002, 2003, papers II 
and III, respectively) we use the density fields of SDSS and 
LCRS galaxies to derive catalogues of clusters and super- 
clusters and to compare samples of density-field defined 
clusters and superclusters with clusters and superclusters 
found with conventional methods. 

In the next section we describe the Early Data Release 
of the SDSS, and the LCRS samples of galaxies used. In 
section 3 we derive the luminosity function for the SDSS 
and LCRS galaxies using distances found for a cosmolog- 
ical model with dark matter and energy. In section 4 we 
calculate the density helds using SDSS and LCRS galaxy 
samples. We analyse our results in section 5, section 6 
brings conclusions. 

2. The data 

2.1. Early Data Release of the SDSS 

The Early Data Release (EDR) of SDSS consists of two 
slices about 2.5 degrees thick and 65 — 90 degrees wide, 
centred on celestial equator, one in the Northern, and the 
other in the Southern Galactic hemisphere (Stoughton et 
al. [2002 ), and contains about 35,000 galaxies with mea- 
sured redshifts. SDSS Catalogue Archive Server was used 
to extract angular positions, Petrosian magnitudes, and 
other available data for all EDR galaxies. From this gen- 
eral sample we obtained the Northern and Southern slice 
samples using redshift interval 1000 < cz < 60000 km s _1 , 
and Petrosian r*— magnitude interval 13.0 < r* < 17.7. 

2.2. LCRS 

The LCRS is an optically selected galaxy redshift sur- 
vey which extends to a redshift of z ~ 0.2 and is com- 
posed of 6 1.5° x 80° slices, 3 in the Northern and 3 in 
the Southern Galactic hemisphere. The survey contains 
26,418 galaxy redshifts, measured via a 50-fibre or 112- 
fibre Multi-Object Spectrograph. The 50-fibre fields have 
nominal apparent magnitude limits of 16.0 < R < 17.3, 
and the 112-fibre fields 15.0 < R < 17.7). A galaxy was 
included in the calculation of the density fields if its flux 
and surface brightness lie within the official survey lim- 
its (see Lin et al. 1996 or Tucker et al. 2000] for details). 
Differences in sampling density and magnitude limits were 
taken into account using statistical weighting of galaxies, 
as will be explained in section 4. 

3. The luminosity function of the SDSS and 
LCRS samples 

Luminosity function is one of the basic characteristics of 
galaxy population which gives us the co-moving number 
density of galaxies per magnitude interval and it is cer- 
tainly necessary to determine it in order to understand 
selection effects in large scale studies. In this section we 
are going to estimate luminosity function for the above 
mentioned SDSS and LCRS samples. 



For the LCRS this analysis has been done for the criti- 
cal matter density case by Lin et al. ( |1996D (L96) but here 
we are going to reconsider it in case f2 m = 0.3, £l\ = 0.7. 
Luminosity function for the SDSS galaxies in all filters 
(it*, g* , r* , i* , z*) has been determined by Blanton et al. 



(2001, hereafter B01) using the commissioning phase data 
which covered approximately 230 square degrees of the 
sky along the Celestial Equator in the region bounded by 
145° < a < 236° and < 1.25°. Here we shall calculate it 
for the SDSS EDR data which covers ~ 400 square de- 
grees of the sky and thus contains approximately twice as 
much information. 

There exists a variety of methods to estimate lumi- 
nosity function e.g. T/Knax method' by Schmidt ( 1968 ), 
pa ramet ric maximum likelihood method by Sand age e t 
al. (|l979| ) (STY79), 'C~ method' by Lynden-Bell ( |l97lD , 
s tepw ise maximum likelihood method by Efstathiou et al. 
( 1988 ) etc. In our case we follow the STY79 approach 
which assumes special parametric type for the differential 
luminosity function c/>(M), namely the one in the form of 
Schechter function (Schechter |l976|) 



(j> (M) dM 
F (M) 



0.4 In (10) cf)*F(M) 

^q-OA(M-M') 



a+1 



exp [-F [M)) dM, 



(1) 



and finds maximum likelihood estimates for the free pa- 
rameters M* and a. Once we have fixed the form for the 
luminosity function we can immediately write down the 
probability that galaxy i, that is observed at the redshift 
Zi, has absolute magnitude Mi as 

P(Mi,Zi 



P l = P{M l | Zi ) 



Mr 



(2) 



/ <p (M) dM 

M mla ( Z< ) 



where the third equality follows from the usual and rather 
crude assumption that luminosity function is independent 
of the spatial location. Here M max (zi) and M m i n (zt) de- 
termine the observational window in absolute magnitude 
at redshift Zi that corresponds to the surveys limiting ap- 
parent magnitudes 77i max and m mul , so 



M n 

M„ 



^ max 



- 25 - 51ogd L (z^ - 

- k(zi)-A(l,b), (3) 

where as usual d^ (z^) denotes the luminosity distance, 
k (z^ is the k-correction and A(l,b) is the absorption 
term. Now we can immediately express the likelihood func- 
tion C for finding N galaxies with absolute magnitudes Mi 
(i = 1 ... N) as 

N 

C = Y[Pi(a,M*) =>ln£(a,M*) = 



JV 

E 



\ncj){Mi 



In 



(p(M)dM 



M mi „( Zi ) 



(4) 
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and the best fitting values for M* and a are found by 
maximising hxL{a,M*). Assuming Gaussianity, the cor- 
responding error ellipses are found as usual 



ln£ = ln£max - 7: Ax 2 , 



(5) 



where A\ 2 is the critical value for the desired confidence 
level for the x 2 distribution with two degrees of freedom 
(Press et al. 1992| ). Let us now discuss some pros and cons 
for the STY79 method. 
Pros: 

• as maximum likelihood method it gives us optimal es- 
timates (i.e. the ones with the minimum variance) for 
the parameters M* and a, 

• STY79 method is insensitive to galaxy clustering and 
density evolution 

... and cons: 

• the goodness of fit of the proposed parametric lumi- 
nosity function (i.e. Schechter function) cannot be as- 
sessed]^, 

• the normalisation <jj* has to be determined using other 
methods. 

The first of these two problems is not really very seri- 
ous since it is well known that the Schechter function has 
almost always served as a very good approximation and 
in our case it is important to mention that earlier stud- 
ies of the SDSS and LCRS data have certainly demon- 
strated this fact (B01, L96). Since our main task here is 
to estimate the selection function, which certainly doesn't 
depend on the absolute normalisation of the luminosity 
function, we are going to determine <j>* and its uncertainty 
only for the SDSS data using rather simple arguments. 
Namely <fi* is obtained so as to produce the observed num- 
ber counts and its error is estimated using re-sampling 
techniques (to be more precise, " Jackknife" method). Here 
we must also correct for the fact that some galaxies are 
missing due to lack of fibres in dense regions, due to spec- 
troscopic failures as well as fibre collisions. These effects 
together give an average sampling rate of 92% (B01). 

In order to calculate the luminosity distance, di,(z) = 
r(z)(l + z) (for coordinate distance r(z) see equation (||)) 
as well as to estimate k-corrections, we have to fix the 
cosmological model. As already mentioned, we take spa- 
tially flat model with O m = 0.3 and Oa = 0.7. To estimate 
k-corrections 



k(z) = 2.51og(l + z) + 
+2.5 log 



J F(X)R (A) XdX 




lF[ T ^)R(X)Xd\ 



(G) 



1 To overcome this problem Efstathiou et al. ( 1988| ) pro- 
posed stepwise maximum likelihood method. 



we use spectral templates F (A) obtained from Shimasaku 
( 1998 ) and SDSS filter response curves R(X) (which in- 
clude atmospheric and CCD responses) as compiled by 
Strauss and Gunn^j. For these template spectra we first 
calculate g* — r* colours for a range of rcdshifts. Now, for 
each observed galaxy, we interpolate on this grid to find 
the intrinsic g* — r* colours (according to which we clas- 
sify our objects) that correspond to the observed colours 
and redshifts. For the same template spectra we calculate 
k-corrections for different redshifts, and now taking into 
account the obtained intrinsic g* — r* colours, we can in- 
terpolate to find the desired k-corrections for each galaxy. 



Finally, we used the Schlegcl et al. (1998) maps to obtain 
the reddening-corrected absolute magnitudes. 

We are going to estimate the luminosity function of 
the SDSS galaxies only for the r* band (i.e. the one with 
the highest sensitivity) using the following constraints for 
the apparent magnitude and recession velocity: 



m r , : 14.5 . . . 17.5, 

v rec : 1000 . . . 60000 km/s. 



Since for the LCRS we have apparent magnitudes only 
in the "hybrid" Kron-Cousins R filter i.e. there is not any 
colour information, we are going to approximate (as was 
also done by L96) spectral energy distributions by power 
law F (y) oc u~ a , assuming average spectral index a ~ 2. 
The filter is "hybrid" since the LCRS photometry was 
obtained through a Gunn r filter, although the calibration 
was done relative to standard stars in the Kron-Cousins 
R-band (for more information on LCRS see e.g. Tucker 
1994). In this case equation M) can be trivially integrated 



to give us 



k(z) = 2.5(a-l)log(l + z) 



(a~2) 



2.51og(l + z). 



(7) 



To estimate Galactic extinction effects for the LCRS mag- 
nitude s we use directly dust maps by Schlegel et al. 
( 1998 ) assuming that these "hybrid" magnitudes differ 
only negligibly from "true" Kron-Cousins R magnitudes 
(see e.g Tucker 1994 ) in which case absorption Ar kc is 
related to the colour excess E (B — V) through a rela- 
tion A Rkc = 2.673E(B - V) . Magnitude limits for the 
LCRS were quoted in the previous section and for the re- 
cession velocity we take an interval 1000 . . .45000 km/s. 
Here, as the photometry is not so precise as in the case 
of the SDSS, we take into account that the observed lu- 
minosity function is a convolution of the 'true' luminosity 
function with the magnitude errors which we assume to 
be Gaussian distributed with the variance a = 0.1 mag- 
nitude (L96). Since in the LCRS first 20% of data was 
obtained using 50-fibre spectrograph and the rest using 
11 2- fibre system, one must take into account field-to-field 
variations in the sampling ratio. Also we correct for the 
surface brightness selection effects and for the apparent 
magnitude and surface brightness incompletenesses. For 
more detailed information on these issues we again refer 
to L96. 
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Fig. 1. la and 2a error ellipses for the SDSS Northern, 
Southern and total samples. 



Fig. 3. Results of the T/Kxax method' for the total SDSS 
sample with the corresponding Poissonian error-bars. The 
best fitting Schechter function obtained via STY79 maxi- 
mum likelihood method is also given by a dashed line. 




M 



Fig. 2. la and 2a error ellipses for the LCRS N50, S50, 
N112, S112 and total samples. 



The results for the STY79 maximum likelihood 
method are given in Tables [l] and | and in Figures and 
^ for the SDSS and LCRS samples, respectively. Slices 
that are located in the Northern and Southern Galactic 
hemispheres are denoted by N and S, and for the LCRS 
the numbers 50 and 112 denote the number of fibres used. 



For the LCRS this kind of division is the same as used 
by L96 and so their Figure 4 (which was obtained as- 
suming critical matter density) can be compared directly 
to Figure || (which assumes f2 m = 0.3, Qa = 0.7). As 
already mentioned, luminosity function normalisation <jf 
was calculated only for the SDSS data, giving us the value 
4>* = (151 ± 6) • 10" 4 h 3 Mpc" 3 . If the analysis was per- 
formed separately for the Northern (N) and Southern (S) 
slices, then the values 0.0140 h 3 Mpc -3 and 0.0158 h 3 
Mpc" 3 were obtained for N and S, respectively. If we as- 
sumed the best values for M* and a obtained for the whole 
sample (i.e. M* = —20.71, a = —1.12) then normalisa- 
tions for N and S were 0.0152 h 3 Mpc" 3 and 0.0149 h 3 
Mpc -3 , respectively. Our results for a and <j>* for the SDSS 
N slice (which covers the commissioning phase data anal- 
ysed by B01) agree well with the values quoted by B01, al- 
though for the Schechter parameter M* our value —20.76 
is somewhat lower, agreeing only marginally with their 
result. 

In Figure |^ we compare the best fitting Schechter 
function obtained via STY79 method with the results of 
the simple T/Vm ax method' along with the corresponding 
Poissonian error-bars. This serves just as an independent 
cross-check of our results, although it is of course well 
known that '1/V^ lax method' gives unbiased results only 



See jittp : // archive .stsci.edu/ sdss/documcnts / response. dat for a homogeneous distribution (Felten 1976). 
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Table 1. Best fitting M* and a for the SDSS samples 

Sample Al* — 5 log h a 

N - (20.76 ± 0.04) -(1.19 ±0.05) 

S - (20.67 ±0.05) -(1.07 ±0.05) 

TOTAL -(20.71 ±0.03) -(1.12 ±0.03) 

Table 2. Best fitting M* and a for the LCRS samples 



Sample 


Al* — 5 log h 


a 


N50 


-(20.33 ±0.12) 


-(0.40 ±0.18) 


S50 


-(20.64 ±0.18) 


- (0.74 ± 0.21) 


N112 


- (20.40 ± 0.05) 


- (0.76 ± 0.07) 


S112 


- (20.40 ± 0.05) 


- (0.70 ± 0.07) 


NS112 


- (20.38 ± 0.04) 


- (0.70 ± 0.04) 


TOTAL 


- (20.40 ± 0.03) 


- (0.69 ± 0.04) 



4. Density field 

We have extracted from the SDSS and LCRS catalogues 
subsamples of galaxies with redshifts z < 0.2. The calcula- 
tion of the density fields consists of three steps: 1) calcula- 
tion of the distance, absolute magnitude, and weight factor 
for each galaxy of the sample; 2) calculation of rectangu- 
lar coordinates of galaxies and rotation of coordinate axes 
in order to minimise projection effects; and 3) smooth- 
ing of the density field using an appropriate kernel and 
smoothing length. 

Observed redshifts were first corrected for the motion 



relative to the CMB dipole (Lineweaver et al. |1996| ). Co- 
moving dista nces o f galaxies were calculated as follows 
(e.g. Peacock ^999| ): 



r(z) = 



i + a, 



(l±.x) 3 -l]} 



dx. 



(8) 



As already mentioned, we used a cosmological model with 
density parameters: matter density Q. m = 0.3, dark energy 
density J1a = 0.7, total density £7o = f2 OT + f^A = 1-0, all 
in units of the critical cosmological density. With these 
parameters the limiting redshift zn m = 0.2 corresponds 
to co-moving distance ru m = 571 h~ x Mpc. Because both 
surveys cover relatively thin slices on sky, we shall cal- 
culate in the following 2-dimensional density fields, pro- 
jecting all galaxies to a plane through the slice. A 2- 
dimcnsional density field of SDSS EDR was calculated also 



by Hoyle et al. (2002), who compared the geometry of the 
large-scale matter distribution with ACDM simulations. 

Our redshift surveys cover a fixed interval in appar- 
ent magnitude which corresponds to a certain range in 
luminosity, L\ and L2, depending on the distance of the 
galaxy. This range, and the observed luminosity of the 
galaxy, L b s , were found as described in the previous sec- 
tion. 

We regard every galaxy as a visible member of a den- 
sity enhancement (group or cluster) within the magnitude 
window of the group. Further we suppose that the lumi- 
nosity function, calculated for the whole sample, can be 



applied also for individual groups and clusters. We shall 
discuss in Paper II problems associated with the calcu- 
lation of the density field and properties of clusters and 
superclusters, found from density field. Based on these as- 
sumptions we calculate the number and luminosity density 
fields. In case of the number density the weighting factor 
is proportional to the inverse of the selection function, 



1 



w 



(9) 



where n oxp (/, D) is the expected number of galaxies in the 
whole luminosity range 

J °°<f>(L)dL 



cxp 



(f,D)=F 



J£<KL)dL 



(10) 



Here F is the field-to-field sampling fraction, <fi(L) is the 
luminosity function, and luminosities L± and L2 corre- 
spond to the observational window of apparent magni- 
tudes at the distance of the galaxy. The fraction F takes 
into account the difference between 50-fibre and 112-fibre 
data and other effects, for details see L96. We assume 
that galaxy luminosities are distributed according to the 
Schechter function (Schechter |1976| ): 



4>(L) oc {L/L*) a exp(-L/L*)d(L/L*), 



(11) 



where a and L* are parameters. The values of these pa- 
rameters used in our calculations are given in Table ^ for 
each of the slices. For the LCRS —6° N slice we have taken 
different values for a and L* since this slice was covered 
entirely with 50-fibre measurements. 

In the case of luminous density every galaxy represents 
a luminosity 



Lfot = L i, s Wl, 



(12) 



where L obs = L 1O°- 4x(Mq - m) is the luminosity of the 
visible galaxy of absolute magnitude M, M & is the ab- 
solute magnitude of the Sun in an appropriate filter. In 
the SDSS r* filter M© = 4.62 (B01) and in the Kron - 
Cousins R-band M© = 4.47 (Binney & Merrifield |1998| ). 
The weight (the ratio of the expected total luminosity to 
the expected luminosity in the visibility window) is 



W L = 



Jp 00 L<j)(L)dL 
J^L^L)dL 



(13) 



We calculated expected luminosities by numerical inte- 
gration using a finite total luminosity interval, low limit 
Lq corresponds to absolute magnitude Mq = —13.0, and 
upper limit Ln rn to magnitude Mu m = —24.5. In the 
case of SDSS we adopt the apparent magnitude window 
mi = 17.70, TO2 = 13.0 in r*-band; in the case of LCRS 
the apparent magnitude window is different for various 
fields and is taken into account using corresponding ta- 
bles. 

This weighting procedure was used also by Tucker et 



al. (200C) in the calculation of total luminosities of groups 
of galaxies. A different procedure was applied by Moore, 
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Table 3. Data on SDSS and LCRS slices 



Sample 


S 


(R.A.) 


AS 


AR.A. 


a 


M* 


Ngal 


SDSS 


0°N 


190.25° 


2.5° 


90.5° 


-1.12 


-20.71 


15221 


SDSS 


0°S 


23.25° 


2.5° 


65.5° 


-1.12 


-20.71 


11864 


LCRS 


-3°N 


191.4° 


1.5° 


81.0° 


-0.69 


-20.40 


3726 


LCRS 


-6°N 


189.8° 


1.5° 


77.9° 


-0.40 


-20.33 


2132 


LCRS 


-12°N 


191.4° 


1.5° 


81.1° 


-0.69 


-20.40 


3961 


LCRS 


-39°S 


12.1° 


1.5° 


113.8° 


-0.69 


-20.40 


3390 


LCRS 


-42°S 


12.2° 


1.5° 


112.5° 


-0.69 


-20.40 


3610 


LCRS 


-45° S 


12.3° 


1.5° 


114.1° 


-0.69 


-20.40 


3289 



Frenk & White ( |1993| ) (adding the expected luminosity 
of faint galaxies outside the observational window). We 
assume that density and luminosity distributions arc in- 
dependent, which leads to a multiplicative correction. 

To find the number and luminosity density fields we 
calculated for all galaxies rectangular equatorial coordi- 
nates as follows: 



of the particle is distributed between neighbouring cells 
using the Gaussian smoothing 



x = r(z c ) cos <5 cos a, 
y = r(z c ) cos <5 sin a, 
z = r(z c ) sin<5, 



(14) 



here z c is the redshift, corrected for the motion relative 
to CMB dipole. Next we rotated this coordinate system 
around the z-axis to obtain a situation where the new 
y-axis were oriented toward the average right ascension 
of each slice (see Table ||) and finally we rotated this 
new system around its x-axis so as to force the average 
z-coordinate to be zero in order to minimise projection 
effects. 

In Table || we also present some of the characteristic 
parameters for each of the studied slice as the average 
declination {5} and right ascension {R.A.), the number of 
galaxies, N ga i, included in the calculations of the density 
fields etc. 

The EDR of SDSS and LCRS consists of essen- 
tially 2-dimensional slices of about 80° wide and about 
450 h,- 1 Mpc depth. The width of slices is 1.5 - 2.5°, thus 
in space slices form thin conical (wedge-like) volumes; 
the thickness of the cones at a characteristic distance 
300 hr 1 Mpc from the observer is only 8 — 12 h^ 1 Mpc. 
Due to the thin shape of slices we shall calculate only a 
2-dimensional density fields. 

To smooth the density field two approaches are often 
used. The first one uses top-hat smoothing as customary 
in N-body calculations. In this case the density in grid 
corners is found by linear interpolation where the 'mass' 
of the particle is divided between grid corners according to 
the distance from the respective corner. Variable smooth- 
ing length can be obtained by changing the size of the grid 
cell. A better and smoother density field can be obtained 
when we use Gaussian smoothing. In this case the 'mass' 



W jg = (27R7, 



2 ) 

.Sill J 



exp 



2a s 2 m 



(15) 



here and x g are positions of the galaxy and the grid 
cell, respectively, and a sm is the smoothing length. 

To calculate the density fields, we have formed a grid 
of cell size 1 hr 1 Mpc. On farther part of the survey only 
very bright galaxies can be observed, the weight factor be- 
comes large and the fields are noisy. To decrease this effect 
we have calculated the density fields for the LCRS only 
for co-moving distance up to ru m < 450 hr 1 Mpc. Very 
close to the observer the density fields are also noisy. Here 
slices become very thin in real space, and only galaxies 
of low luminosity fall into the window of apparent magni- 
tudes observed. For nearby region the density fields were 
calculated but density enhancements were not considered 
as clusters if r < 50 hr 1 Mpc. 

In calculations we have used several smoothing lengths. 
First, we tried a smoothing length <r sm = 2 h^ 1 Mpc. This 
smoothing yields a very smooth field where density en- 
hancements (clusters and groups of galaxies) can be easily 
identified. However, the size of these enhancements clearly 
exceeds the size of real clusters of galaxies. Thus we have 
used a smaller smoothing length, <r sm = 0.8 h^ 1 Mpc, 
to calculate a density field which can be considered as 
an approximation to the true density field of dark matter 
associated with luminous galaxies. 

To identify superclusters of galaxies we have used a 
larger smoothing length, <7 sm — 10 /i _1 Mpc. Experience 
with numerical simulations has shown that this smooth- 
ing length is suitable to se lect su percluster-size density en- 
hancements (Frisch et al. 



1995). To take into account the 



thickness of the slice, the smoothed density fields were di- 
vided by the volume of the column at the location of the 
particular cell in the real 3-D space, covered by observa- 
tions of the slice. In this way the density fields are reduced 
to a plane parallel sheet of constant thickness. 

The high-resolution number and luminosity density 
fields for the SDSS EDR Southern slice are shown in 
Figure g. The density fiel ds of all EDR SDSS slices are 
presented in Einasto et al. (2002), the fields of LCRS slices 
are presented in Einasto et al. (2003). 
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Fig. 4. Density fields of the SDSS EDR Southern slice, smoothed with a = 0.8 h 1 Mpc dispersion. Densities are 
expressed in units of the mean density. Upper panel shows the number density field, lower panel the luminosity density 
fipld 
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5. Analysis 

So far we have determined the number and luminosity den 
sity fields for the SDSS EDR and LCRS samples. In orde 
to estimate the selection effects we have calculated the lu 
minosity function assuming Schechter type of parametri 
sation. 

To be sure that we are not grossly under- or overesti 
mating the statistical weights we are applying to visibl 
objects (in order to take account their invisible compan 
ions), we perform some simple calculations. Namely, wi 
calculate the number and luminosity density averaged ove 
thin shells as a function of the coordinate distance. Wha 
concerns the number density, then we expect it to fall h 
the distant part of the samples since there we miss some o 
the faintest systems totally. The same need not be true for 
the luminosity density since even the modest amount of 
luminosity evolution could compensate the decrease in the 
number density. The results for the number and luminos- 
ity density for the SDSS EDR data are given in Figure pi. 
We see that indeed the number density starts to decline at 
the largest distances, as expected, indicating that we are 
at the right track. Contrary to the number density the lu- 
minosity density is almost constant for the Northern slice, 
but shows a slightly increasing trend in the Southern slice 
on large distances. This is also evident from the Figure 
H where densities are colour coded. For the LCRS slices 
we get essentially the same results and for the brevity we 
are not going to present them here. All this suggests that 
our statistical weighting scheme looks quite reasonable for 
such a preliminary study. 

The second thing we want to mention is the fact 
that we have assumed that the same Schechter param- 
eters M* and a apply for the whole sample. In order to 
study the possible evolution with redshift z we have di- 
vided our SDSS EDR sample into two parts: low z part 
(v rec — 1000... 30000 km/s) and high z part (v rec — 
30000... 57000 km/s) and performed the STY79 maxi- 
mum likelihood analysis for both of these separately. The 
results that are given in Figure ^ probably suggest that 
M* and a are evolving in time. From this figure we obtain 
a simple linear approximation 

M* = -4.2z- 20.37, 
a = -6.7z-0.63, (16) 

where we have taken into account the fact that the aver- 
age z for the low z sample is 0.067 and for the high z one 
0.136. If we assume this kind of M* and a dependences 
on z for the whole sample and calculate the number and 
luminosity densities as above then we get pretty bad re- 
sults. So certainly this kind of big linear variation with z 
makes things worse. 

Next we divide the SDSS EDR Southern sample 
into three parts: v rec — 1000... 15000 km/s, v rec = 
15000 . . . 30000 km/s and v rec = 30000 . . . 60000 km/s. In 
this case we get non-monotonic dependences for M* and 
a on z (see Figure]^). We also notice that variations in M* 
and a are correlated with the density of the surroundings: 




Fig. 5. Upper and lower panels show respectively the 
number and luminosity density as a function of coordi- 
nate distance for the SDSS Northern and Southern slices 
using the best fitting Schechter parameters for the total 
sample as given in Table 0. Densities were calculated using 
the high-resolution field. 
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Fig. 6. lcr and 2a error ellipses for the SDSS data divided 
into low z and high z subsamples. 
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Fig. 7. la and 2a error ellipses for the SDSS EDR 
Southern slice data divided into three z intervals. 
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Fig. 8. la and 2a error ellipses for the SDSS EDR 
Southern slice data divided into three global density S 
intervals. 



due to the effects of the nearby dense regions the value 
for M* is smaller than on average; as we move further 
away M* gets bigger, which is caused by the large void 
at distances ~ 200 . . . 350 h~ x Mpc (see Figure |). In dis- 
tant parts, as the density rises again, M* gets accordingly 
smaller. This suggests that environmental effects are more 



Table 4. Luminosity function parameters in different en- 
vironment 



Sample 


M* — 5 log h 


a 


< So < 1.0 


-20.55 ±0.07 


-1.14 ±0.06 


1.0 < So < 2.0 


-20.67 ±0.07 


-1.09 ±0.05 


2.0 < So 


-20.81 ±0.07 


-1.13 ±0.05 
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Fig. 9. Luminosities of Southern slice galaxies as a func- 
tion of the global relative density S smoothed with 
2 h^ 1 Mpc dispersion. Luminosities are expressed in units 
of 10 10 Solar luminosities, global density g in units of the 
mean density. The Northern slice has similar dependence 
on the environmental density. 



important than evolutionary ones (which we expect to be 
monotonic in z), and so the usual assumption in determin- 
ing the luminosity function (f> (M) that it is independent 
of the spatial position cannot be strictly correct. 

To study this aspect in more detail we calculated the 
dependence of luminosity function parameters on the den- 
sity of the environment. We calculated for every galaxy 
the global density So (with 10 /i _1 Mpc smoothing) and 
used this value as environmental parameter. Galaxies were 
divided into three types according to density (expressed 
in units of the mean density). Parameters of the lumi- 
nosity function are given in Table ^ and Figure ||. We 
see that parameter M* changes considerably with den- 
sity: in high-density regions its value is more negative, 
i.e. galaxies are brighter. A similar tendency has been re- 



cently observed in other surveys (Bromley et al. (1998) 



Beisbart k Kerschcr (2000)| , |Norberg et al. (2001)| ). 

The density field with 10 ft. -1 Mpc smoothing charac- 
terizes the density on supercluster scales. The dependence 
of galaxy luminosity on density has been found also on 
smaller scales. In Figure M we plot luminosities of galaxies 
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of the Southern SDSS EDR slice as a function of the den- 
sity smoothed on 2 hr 1 Mpc scale. This smoothing is sen- 
sitive to density changes within clusters of galaxies and in 
galaxy filaments. We see a strong dependence of the upper 
end of luminosities: brightest galaxies in low-density re- 
gions are about 5 times less luminous than in high-density 
regions. This difference in luminosity corresponds to a dif- 
ference 1.7 magnitudes, much more than expected from 
the analysis based on larger smoothing length. We plan to 
address all these issues in more detail in the forthcoming 
papers. Physical interpretation of these phenomena needs 
theoretical studies. 

We repeat that all the density fields (see Figure ||) 
compiled in this paper refer to the redshift space. 



6. Conclusions 

In this paper we have calculated the galaxy luminosity 
function for the SDSS EDR and LCRS samples and used 
it to construct the number and luminosity density fields 
(smoothed on 0.8 h~ l Mpc scale) assuming flat underlying 
cosmologies with Cl m = 0.3 and J7a = 0.7. The analysis 
presented here is rather simple and serves as a first step 
in the study of the distribution of galaxies in SDSS and 
LCRS samples. The principal conclusion from our study 
is: parameters of the galaxy luminosity function depend 
on the distance from the observer, density of the environ- 
ment, they are different for the Northern and Southern 
slice. The largest effect is the dependence on the density of 
the environment: in high-density regions brightest galaxies 
are more luminous than in low-density regions by a factor 
up to 5 (1.7 magnitudes). Some of these effects suggest 
that it is not yet possible to find an universal set of pa- 
rameters of the luminosity function valid for a fair sample 
of the Universe. In other words, presently available sam- 
ples are still too small to be considered as candidates for 
the fair sample. 
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